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Abstract 

We re-examine the effort to constrain the time- variability of the coupling 
Qh! constants of the fundamental interactions by studying the anomalous iso- 

r-| ! topic abundance of Sm observed at the remnants of the natural reactors 

which were in operation at Oklo about 2 billion years ago, in terms of a 
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possible deviation of the resonance energy from the value observed today 
We rely on new samples that were carefully collected to minimize natural 
contamination and also on a careful temperature estimate of the reactors. 
We obtain the upper bound (—0.2 ±0.8) x 10~ 17 y _1 on the fractional rate 
of change of the electromagnetic as well as the strong interaction cou- 
pling constants. Our result basically agrees with and even suggests some 
improvement of the result due recently to Damour and Dyson. Strictly 
speaking, however, we find another range of the resonance energy shift 
indicating a nonzero time variation of the constants. We find a rather 
strong but still tentative indication that this non-null range can be ruled 
out by including the Gd data, for which it is essential to take the effect 
of contamination into account. 



*A preliminary report is found in Ref. [1]. 



1 Introduction 



Few aspects of physics can be understood without recourse to fundamental "constants," 
such as the speed of light c, Planck's constant h, Newton's constant G, the fine-structure 
constant a, and its QCD counterpart «qcd, or the pion-nucleon coupling constant squared 
a s = gl N /4:TT in most of realistic analyses. To this list one may also add the masses of 
certain fundamental particles, like the quarks and the leptons, or some of the vacuum 
expectation values, and perhaps the cosmological constant A. The traditional assumption 
that these parameters of physics theories are constant was challenged for the first time 
by Dirac, who suggested that G might vary as t -1 with t the cosmic time [2]. This paper 
inspired to many probes of the possible time variability of various fundamental constants. 

However, no positive evidence for any time-dependence has been discovered so far, nor 
does Dirac's original argument seem fully convincing. Nevertheless, the very idea that at 
least some of the fundamental constants may not be truly constant still attracts serious 
attention partly because this appears to be a rather general outcome of the recent efforts 
toward unifying particle physics and gravity. 

Quite remarkable is the finding that the upper bounds obtained so far for a and a s 
are many orders of magnitude below the value ~ 10~ 10 y -1 expected naturally in terms 
of the present age of the universe to ~ 10 10 y [3-9]. The same tendency is indicated also 
for G [10]. This is expected to provide an important clue to the nature of a model of 
unification. 

The most stringent constraint ever reported is due to Shlyakhter [6,7] who exploited 
the exceptionally sensitive dependence of the cross section of the reaction 

n + 149 Sm -> 150 Sm + 7, (1) 

on the energy of a resonance lying as low as E r = 97.3 meV above threshold, which 
corresponds to the temperature 856° C and is much smaller than any of the mass scales 
of the strong or the electromagnetic interaction. 

It was observed that the anomalous abundance of 149 Sm observed at Oklo [11,12] can be 
understood rather well in terms of the reaction (|1]) with the fundamental-constant values 
observed today. An analysis of the uncertainties in the Oklo observations provides an 
upper bound on the possible deviation of the neutron-capture cross section a 149 associated 
with the reaction (0) occurring 2 billions years ago from its present value. Shlyakhter 
concluded that | Aaug/oug | <,10%, which implies \AE r \ <,20 meV. Furthermore he derived 
that |Aa s /a s | <5 x 10~ 10 y _1 , or |d s /a s | <,2.5 x 10~ 19 y" 1 , which represent upper bounds 
on the time-variability of the fundamental constant that are several orders of magnitude 
more stringent than any other estimate of these bounds [3-5,8,9]. 

However, this substantial improvement of the upper-bound estimate does not seem 
fully appreciated, because complete details of the analysis have never been published. 
Another reason for the limited acceptance may be that the derivation steps are much less 
direct than in estimates based on quasi-stellar-object spectra [8] or on clock standards [9]. 

More recently Damour and Dyson [13] repeated the analysis on the same process ([J) 
in greater detail, reaching a somewhat more conservative estimate, — 120meV < AE V < 
90meV. They also showed that the result can be interpreted as an upper bound \a/a\ <6 x 
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10 y , through the time variation of the Coulomb energy which can be analyzed on a 
much more solid theoretical basis than for the strong- interaction coupling constant a s . 

Having learned, however, about the availability of the data taken even more recently, 
we decided to pursue the subject further. The data are samples obtained from deep 
underground with a geologist's expertise so as to limit potential outside contamination as 
much as possible. We also try to use the data from other isotopes to confirm the result 
obtained from Sm. Possible check of the effect of outside contamination is also planned. 

We find two ranges for the allowed values of AE T , one that covers the null result and 
another apparently indicating a nonzero effect. We try to see if we can exclude one of them 
by the aid of the Gd data which is admittedly subject to more uncertainties than the Sm 
data, due to much smaller abundance because of the more significant absorption effect. 
However, by taking advantage of the presence of the two low-lying isotope resonances, we 
are able to suggest that the range around AE T = is favored if we accept certain plau- 
sibility arguments. A narrower constraint on the reactor temperatures obtained recently 
provides an added support to this conclusion. 

We also show that the theoretical analysis designed originally to be applied to con- 
strain the variability of the electromagnetic coupling constant [13] can be applied to the 
strong interaction as well, thanks to the exceptionally small value of E T , requiring a near 
cancellation between the effects of the two interactions. 

We introduce basic equations that are used in our analyses: in Section 2 those related 
to neutron absorption and in Section 3 those related to nuclide transformations along 
isotope chains. New measurements on five relatively recent samples are then presented 
and analyzed in Section 4, providing a better upper bound on AE r . Section 5 discusses the 
relevance of the results to the possible time-variation of the electromagnetic-interaction 
and the strong-interaction coupling constants. Section 6 gives discussions summary of the 
results obtained in our investigation. Appendix A provides some details of the numerical 
calculation of the thermal average of the resonance cross sections. In Appendix B we 
show some details of the fluence determination adopted in this investigation. 



2 Neutron absorption cross sections 

The results obtained in this investigation are based on an analysis of the neutron-absorption 
cross section in the reaction ([[]) and the corresponding reactions on Gd. The reactions 
are dominated by a resonance, which we describe in the standard Breit-Wigner form, 

_ goTfh 2 r n r 7 

av ~ 2mE(E-E T r + rlj4> [ } 

where g = (2 J+l)/[(2s+l)(2/+l)] is the statistical factor, s the spin of the neutron, / the 
spin of the target, J the spin of the compound nucleus, and r tot = T n + r 7 the total width 
in terms of the neutron and the photon widths, respectively. Also, T n = (2k/K)(D/ir), 
where k 2 = 2mE/h 2 , K the wave number inside the target nucleus, and D the level 
spacing of the compound levels. The values of these parameters are given in Table [1] for 
the lowest-lying resonance in the reactions on 149 Sm, 155 Gd and 157 Gd. 
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Channel 


149g m 155Q d 157Q d 


Resonance energy E t (meV) 
Neutron width T n (meV) 
Gamma width T 1 (meV) 
Statistical factor go 


97.3 26.8 31.4 
0.533 0.104 0.470 

60.5 108 106 
9/16 5/8 5/8 



Table 1: Parameters of the lowest-lying resonances used in the present analysis. The tabulated 
values of E T are the present resonance-position values E t q [14]. The listed values of T n are 
estimated at E T . 

We now derive how the isotopic abundances observed today depend on E T , since we 
need to consider the possibility that E T were different from its present value 2 billion years 
ago. Other parameters may be assumed constant, since only E T will affect the result in 
any significant manner. Note that we would normally expect abundance changes not 
much larger than Act/ a, unless they are amplified by the resonance mechanism and by 
the extremely small value of E T relative to any of the mass scales of nuclear physics. 

We average (Q) with respect to the thermal neutron flux, 

<$>{E,T) =v{E)p{E,T), with v = ^2E/m, (3) 

where 

p(E, T) = 2T,- 1 ' 2 T-^ 2 e~ E / T ^E (4) 

is the normalized Maxwell-Boltzmann distribution. 

Normally the thermally averaged cross section is defined by 

o(T) = [ViT)}- 1 J a(E)$(E, T)dE, (5) 

where 

V{T) = J $(E,T)dE. (6) 

However, in most of the analyses of the Oklo phenomenon, it is customary [15] to 
replace the denominator T>(T) by 

V = v Jp(E,T)dE = v , (7) 

for a convenient choice of the velocity Vq = 2200 m/s corresponding to the temperature 
T = {m/2k)v 2 = 20.4°C, hence defining an effective cross section 

a(T)=V- 1 J a t (EME,T)dE=^^a(T). (8) 

The advantage of this non-standard definition is that, unlike T>(T), D is a constant, 
independent of T. We note that the neutron flux at Oklo 2 billion years ago is determined 
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Figure 1: Calculated thermally averaged cross sections cii49 for n + 149 Sm — > 150 Sm + 7 as 
functions of the resonance position change AE r and the temperature T. 

based on the observed abundances of 143,145 Nd and 147 Sm, as will be shown in more detail 
in Appendix B. The corresponding neutron-absorption cross sections are known to obey a 
1/v behavior to a good approximation, implying that the effective cross sections of these 
processes are nearly T- independent, hence simplifying the practical analysis considerably 
In the following analysis, the cross sections always occur multiplied by the thermally 
averaged flux </>(T). It is obviously convenient to introduce an effective flux in such a 
way that the product with a remains the same as the original product: 

<j(f) = acf). (9) 

Using the last expression of (|[), we find 

# = ^ do) 

We evaluate the integral in (pD by numerical integration. Details of the calculation 
are found in Appendix A. The results for reactions on 149 Sm, 155 Gd, and 157 Gd are shown 
in Figs. respectively. The figures illustrate the strong dependence of the thermally- 
averaged cross section on the resonance-position change AE r for various reactor temper- 
atures T. 



5 




-200-150-100 -50 50 100 150 200 
Resonance Position Change AE r (meV) 

Figure 2: Calculated thermally averaged cross sections 17155 f° r n + 155 Gd -^ 156 Gd + 7 as 
functions of the resonance position change AE T and the temperature T. 

We have calculated numerically the effect of thermal fluctuations of the atoms in 
the neutron absorber on the calculated cross sections. The magnitude of this so-called 
Doppler effect on the calculated <x(T) is less than 1% even at the highest temperature 
investigated, 1000°C, and much smaller for lower temperatures. Therefore this effect can 
be safely ignored, and the formalism below is for simplicity developed without including 
Doppler broadening. 

3 Relations for isotopic compositions 

The isotope ratios present in Oklo today are superpositions of contributions from fission, 
subsequent neutron capture, and of the original natural abundances. The interpretation 
of the observed abundances is simplified when in special cases some processes do not 
contribute. In a few of those cases, where the various contributions to the now observed 
isotope ratios are most easily disentangled, we may draw stringent conclusions about the 
variation in time of a few low-lying neutron-resonance positions. We consider here the 
cases of neutron capture on 149 Sm, 155 Gd, and 157 Gd. 
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Figure 3: Calculated thermally averaged cross sections 17157 f° r n + 157 Gd -^ 158 Gd + 7 as 
functions of the resonance position change AE T and the temperature T. 

3.1 Isotopic composition of Sm 

We can make a few observations that simplify the differential equations obtained for the 
time evolution of the number of atoms per unit volume N^(t) of the samarium isotopes 
^Sm near 149 Sm that are connected by nuclear transformations. The spontaneous-fission 
half-lives of 235 U and 238 U exceed 10 16 y. Therefore fission products are generated in 
appreciable amounts only by neutron-induced fission during the period of operation of 
the natural reactors at Oklo. Some isotopes are shielded by stable nuclides and cannot 
be reached in the /3-decay of the unstable nuclei created during the fission process. In 
addition, it is possible to neglect some neutron-capture reactions because of their low 
cross sections, 2.413 barn [16], for example, for capture on 148 Sm. 
Thus, we obtain the following differential equations: 

-a 147 (j)N 147 (t) + N% 35 exp(-o- a 0t)o- f23 50'i / i47, 
6 147 (j)N 147 (t), 

-a U9 (f>N 149 (t) + N$ 35 exp(-a a 0t)a f23 5</>Yi49, (H) 



dN 147 (t) 
dt 

dN 148 (t) 
dt 

dN 149 (t) 
dt 
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where <jf235 is the 235 U fission cross section, a a is the effective neutron-absorption cross 
section of 235 U which will be defined below, Y" 147 and Y li9 are the fission yields of 147 Sm and 
149 Sm, respectively, and the effective flux defined in Eq. (1C). The effective neutron- 
absorption cross section a a is not the bare fission cross section but an effective cross 
section which includes the restitution of 235 U. Restitution refers to the sequence of events 
in which a neutron is captured by 238 U and 235 U is the final decay product created via 
successive /3-decays yielding 239 Pu and a subsequent a decay. The way we determined the 
effective flux cf) and the effective neutron-absorption cross section <j a is given in Appendix 
B where these two quantities and the epi-thermal index are obtained simultaneously by 
solving three coupled equations. The effective neutron-absorption cross section is defined 
in terms of the 235 U total neutron-absorption cross section given by 0235 as 



0" a 



T — Crz)o"235- 



(12) 



The interpretation of this definition is clear when we refer the third of the four equations 
just above equation (|B.12|) . The time t starts at the beginning of the reactor period, with 
4> assumed to be constant in time.[| Obviously 0149 is dominated by the thermally averaged 
cross section of the absorption process (P while cr 147 is the non-resonant absorption cross 
section of 147 Sm. 

The initial number densities iV^(O) are N® 35 for 235 U and 



N A (0) 



jy na t^j>nat 



A 5 



(13) 



for the samarium isotopes, for which the relative fractional natural abundances have 
been observed to be [16]: 



pnat 
- ft 144 



0.0310, R 



nat 
147 



0.1500, R 



nat 
118 



0.1130, and R 



nat 
149 



0.1380. 



(14) 



The value for the stable and inactive 144 Sm is included, because it is needed later, whereas 
the final results will be independent of the overall normalization constant jV nat and the 
capture cross section &147. 

The solution of the system of differential equations flTT| ) is: 

^235^235 ^149 



N 14 9(t) 



exp(-<xi4 9 0t) -exp(-a a 0t) +N mA R% exp(-6" a 0t) 



N M7 (t) + N U8 {t) 



0"149 

^235^ / 147<7f235 



1 - exp(-ajt)] + iV nat (Rl% + R 



148 ) 1 



(15) 



from which it follows that 

AW*i) + AW*i) 



Q"f235fey"l47 



exp(-o- a 0ri y 



+ (R 



nat 1 pnatN 
147 ' rt 148j 



fc0"f235^149 
5"a — Cl49 



exp(-cri490iij - exp(-cr a 



I Dnat 



(16) 



i'lgexp^-o-a^tij 



T Throughout this work, we used this assumption. However, if we use the technique used in Ref. [13], 
it is possible to remove this assumption. Then <j>t should be replaced J <pdt in all our expressions. The 
numerical results given later are unchanged even if we use this generalized expression. 
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at t\, the end of the reactor activity. Since then the samarium isotope concentrations 
N^it) have remained unchanged until today, unless corresponding, naturally-occurring 
isotopes have flowed into the sample regions from the outside. 

The ratio k = N2 35 /N na,t can be determined in the following way. We observe that 
144 Sm is never produced as a fission product and has a negligible capture cross section, 
1.64 barn [16], for neutron capture, that is A^ 144 (t) is constant. We may then calculate 
the ratio iV 1 44(i 1 )/(iV 1 4 7 (t 1 ) + A^" 148 (t 1 )) and compare to the observed ratio: 



R: 



o 



Observed 



#144 (*l) 



L#147(tl) + #148(tl)J CaJculated 



jynat^nat 



#235^147%^ [ 1 - exp(-cr a 0ti) ] + N^(R\% + Rftl] 

u a 



pnat 
-fl-144 



U a. 



+ (R\t\ + R nai 



(17) 



148 ) 



From these equations, we obtain 



pnat 
rl 144 



D / onat i pnat > 
-KO 1-^147 + -^148; 



RnYi 



O 1 147 



0"f235 



[ 1 - exp(-cr a 



which can be estimated from the observed values only. 

We finally solve (|i~6|) for a^g by using observed abundances, and the fluence 4>t\ and 
the effective neutron-absorption cross section <r a defined earlier and discussed in Appendix 
B. The values of fluences and effective cross sections are tabulated in Ref. [17] and using 
them, we obtain the results for five recently analyzed samples shown in Table . 



3.2 Isotopic composition of Gd 

For gadolinium we may define relative fractional natural abundances R™ 1 in complete 
analogy with samarium in Eq. fll~3|). For Gd the relative fractional natural abundances 
for the relevant isotopes are 

R™1 = 0.1480, R™1 = 0.2047, = 0.1565, 

R™1 = 0.2484, and R^ Q = 0.2186. ^ 

Because thermal neutron capture cross section of 160 Gd is 0.796 barn [16], we can use the 
value of the nearly constant Ni 60 (t) observed at each sample location to determine N nBjt . 
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SF84-1469 


SF84-1480 


Sample 
SF84-1485 


SF84-1492 


SD.37 


<f>t x (1/kb) 


0.525 


0.798 


0.622 


0.564 


0.780 


Nxu(tx) (%) 


0.1052 


0.2401 


0.2073 


0.1619 


0.06909 


NxM (%) 


55.34 


53.23 


54.03 


54.81 


52.74 


Nx4s(tx) (%) 


2.796 


3.468 


3.079 


2.890 


4.694 


iV 149 (ti) (%) 


0.5544 


0.2821 


0.4466 


0.4296 


0.3088 


N 2 3 5 (tx)/N 238 (tx) 


0.03181 


0.02665 


0.02971 


0.03047 


0.02435 


(kb) 


85.6 


96.5 


83.8 


99.0 


89.5 



Table 2: Measured isotopic ratios and nuences for Sm for five samples and corresponding, 
calculated values of the cross section 0149. 

Because the neutron-absorption cross sections 0156 and ax58 are 2.188 barn and 2.496 
barn [16], respectively, and because 154 Gd is shielded by the stable 154 Sm the study of the 
possible shifts of the low-lying resonances in 155 Gd and 157 Gd is considerably simplified: we 
may study the transmutations related to the isotope pairs 155 Gd- 156 Gd and 157 Gd- 158 Gd 
separately. The differential equations for the number of atoms N^it) per unit volume for 
the first isotope pair are: 



dN 155 (t) 
dt 

dN 15(i (t) 
dt 



-0"155 0-^155 (*) + ^235 exp(-£ a 0t)a m 50Yi 



155? 



+o 155 0A^i 55 (t) + iV 2 3 5 exp(-a a 0t)of2350yi 



156- 



(20) 



where we ignored the small cross section 0156 = 2.188 barn. We find the solution: 

-^235^/235^155 



AW*) 



N 155 (t) + N 156 (t) 



In analogy with Eq. (K 



Oa - O155 
rnat rynat 



exp(-oi 55 0t) - exp(-o a 0t) 



+iV nat ^exp(-o 155 ^) 

<5"f235 



^235 (^155 + ^156 



1 - exp(-o a 0t) 



we introduce 



+ R 



156 



R, 



o 



Nxm(tx) 



Nx 55 (tx) + N 156 {tx) 



(21) 



(22) 



which can be compared with the observed data expressed in terms of N 160 /N 156 and 
N 155 /N 156 with reference to the stable isotope 160 Gd. In analogy with Eq. ([18]) the ratio 
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k = A^2 35 /A^ nat is here given by 



k 



pnat p ( pnat _j_ pnat\ 
-"■160 ~ K 0{ttl55 + -"156 J 



Ro ?jm (y 155 + r 156 ) [l - exp(-<3- a 0t 1 ) 



Now 6" 155 can be determined by solving iteratively 



N 155 {h 



A 



W 156 (ti) 5' 

where the observed ratio is substituted into the left-hand side, and where 

^5-f235^155 



A = 

B = 



0"a — 0"l55 

fc^l55 O" f 2 35 



exp^-aiss^tij - cxp{-a a (pt 1/ 



°"l55 ~ °a 0"a 
Cf235 



°"155 



1 - exp(-o- a 0ii) 



(7; 



+ exp(-a 155 0t i ; 

1 - exp(-CTi 55 0ii) 



156' 



1 - exp(-<J a 0t 1 ) 1 - exp(-<7i 55 <^l) +B156. 



(23) 



(24) 



(25) 



The transmutation equations related to the isotope pairs 157 Gd and 158 Gd are obtained 
by replacing quantities related to 155 Gd and 156 Gd by quantities related to 157 Gd and 158 Gd, 
respectively, in the above equations. 



4 Analysis of the data 

We use the isotopic data of five samples recently taken from reactor zones (hereafter called 
RZs) 10 and 13 which are located deep underground, 150-250 m below the surface, while 
the other RZs 1-9 are located near the surface. The use of samples from deep-lying RZs 
is important for this study, because it minimizes contamination from in-flow of isotopes 
from the outside after the cessation of reactor operation. Particularly relevant factors 
that influence the contamination process are: 

1. Hexavalent uranium U(VI) species are water-soluble, while tetravalent species are 
insoluble. The valence of uranium is largely determined by oxidation conditions in 
nature. 

2. A uraninite sample located near the surface easily forms water-soluble species be- 
cause of contact with the oxidizing atmosphere. On the other hand, a sample deep 
underground remains insoluble because of a relatively reducing environment. 

3. Uranium minerals under oxidizing conditions are dissolved and recrystalized. It is 
possible that rare-earth elements originally incorporated into uranium minerals are 
released as the minerals are dissolved and precipitated on rock surfaces. 

It is obvious that any geological alteration of the RZs would reduce the quality of the 
isotopic data. Rare-earth-element enriched regions were actually found in layers around 
RZs 1 to 9 [18,19]. For example, the Sm content around GL13 at RZ 9 shows a variation 
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in the range from 21 to 285 ppm, a marked difference from the range 0.43 to 3.32 ppm 
observed in SF84 which is close to RZ 10 used in this study. This geological observation 
suggests that we as a precautionary measure, to keep the contamination to a minimum, 
collect samples from deep-lying RZs. The quality of our new data obtained in this way is, 
for example, demonstrated by a much smaller scatter of the values we obtained for &ug 
compared to what was obtained in Ref. [13] from earlier samples taken from RZs 1-9. 

4.1 Description of the samples 

Because near the surface, rare-earth elements and U often fail to form closed systems, 
they are likely get mixed with outside, normal isotope distributions. This may yield an 
estimate of the fluence that is either too large or too small. We selected our samples 
in order to minimize the effect due to this process. Four samples from SF84 were taken 
from a single drill core which crosscuts RZ 10 located 150 m deeper than the RZs near 
the surface. The SD.37 sample was from a part which was highly depleted in 235 U, in a 
borehole which crosscuts RZ 13. Detailed characterization of the reactor was performed 
through chemical and isotopic analyses of fission products and U taken every few cm from 
the drill core [17]. 

4.2 Sm data 

Table |2| shows measured and calculated results for the five samples labeled in the first 
row. The first four samples were taken from RZ 10, whereas the last sample was taken 
from RZ 13. Measurement errors are given in Ref. [17]. Because they are mostly below 
the level of percent they are negligible in our analysis. More details on this question will 
be reported elsewhere. The last row shows cr 149 calculated as described in Subsection 3.1. 
The result can be represented by 

(Jug = 91 ± 6 kb. (26) 

We then use Fig. [I] to estimate AE T = E T — E r0 , for the assumed temperatures, with the 
results shown in Fig. |4]. 

Note that for the temperature range of interest there are two separated regions for the 
allowed values AE r for a given set of &143 and the lower and the upper temperature limits 
Slower and T upper . In Fig. 1, these regions are bounded by both ends of (p6j), Ti ower and 
^upper, and the left-branch of the curves for one while the right-branch for the other. The 
right-branch range passes through the area that covers AE r = whereas the left-branch 
range extends to the far- negative region of AE r . 

By taking advantage of the strong temperature-dependence of the neutron-capture 
cross section of 176 Lu, most estimates of the reactor temperature have been made based 
on observed isotopic ratios of 176 Lu/ 175 Lu. It has been suspected, however, that the result 
is highly problematic. In one of the careful analyses the temperature range 200 — 360°C 
was obtained [20] , whereas the same technique applied to another set of samples produced 
380°C for one sample with all the six others yielding > 1000°C [17], indicating still un- 
known flaws, such as the lack of more precisely determined temperature-dependence of the 
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Figure 4: The ranges of AE T , represented by horizontal lines, calculated from Fig. [I] corre- 
sponding to the estimated cross section a^g given by (25) are shown for several choices of the 
temperature, namely T = 20, 100, 200, 300, 400, 500°C. 

capture cross section particularly at higher temperatures, possible outside contamination, 
and others, depending on different conditions of reactors. 

In view of this complicated and even confusing situation, it might be worth looking for 
other types of approach. A range of 300-350°C was obtained, for example, by measuring 
the oxygen isotopes of chiorites [21] outside reactor cores, though less direct. This is sup- 
ported by the most recent analysis by Pourcelot and Gauthier-Lafaye who have obtained 
the reactor temperature from the crystallization temperature of chiorites in the reactors 
[22]. Based on this result, and by considering other possible factors, like distance to the 
reactor core and the cooling efficiency, the authors concluded 200-400°C, which we adopt 
here as the most reliable temperature range currently available. Notice that this range is 
consistent with the Lu result in Ref. [20]. 

We then obtain 



AE r = 4 ± 16 meV, for T = 200-400°C, 
for the right-branch range, while 

AE r = -97 ± 8 meV, for T = 200-400 °C, 



(27) 



(2c 



for the left-branch range. The standard deviation has been determined from the lower 
and upper limits of the temperature. 
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c 


SF84-1469 


SF84-1480 


SF84-1485 


SF84-1492 


SD.37 


0.00 


85.6 


96.5 


83.8 


99.0 


89.5 


0.01 


86.4 


100.2 


85.6 


100.7 


90.4 


0.02 


87.1 


104.3 


87.4 


102.4 


91.3 


0.03 


87.8 


108.6 


89.3 


104.2 


92.2 


0.04 


88.6 


113.4 


91.2 


106.0 


93.2 



Table 3: Estimated cross sections CT149 (kb) for five values of the post-reactor contamination £ 
for five samples. 

We also notice that restricting the temperature range further to 300-400 °C, for ex- 
amples, as suggested by Refs. [20,22], does not improve the result ( |27D in any significant 
manner; 

AE T = 5 ± 16 meV, for T = 300-400°C. (29) 

We point out that the shift as large as ~ — lOOmeV suggested in ( f28"D would imply 
a level even below the threshold, hence requiring a more careful analysis. It seems still 
unlikely that it would change substantially the claim that E T at the time of the natural 
reactor was different from the present value beyond many standard deviations. 

The bound thus obtained for AE r , particularly for the right-branch, are rather in- 
sensitive to an increase of T upper as long as T upper <T CTit « 820°C, a critical temperature 
above which the two regions are bridged to each other. In fact for T upper = 700°C, for 
example, we find —83 ± 22meV, and 11 ± lOmeV for the left-branch and the right-branch 
ranges, respectively. The separation between the centers of the two ranges is 94meV, 
which is nearly three times as large compared with the sum 32meV of the half-widths 
of each region. Even if we take the two standard deviation bound for a^g, the critical 
temperature is estimated to be T crit « 730°C, and the ranges for T upper = 700°C are found 
to be —78 ± 30meV and 6 ± 26meV, respectively, which are still well separated, with the 
separation 84meV and the sum of the half-widths 56meV. 

It seems also useful to note that assuming temperatures ,> 1000°C leaves little room 
for AE r in Fig. [I] for (|26|), making higher temperatures less likely from the point of view 
of theoretical analyses. 

As we alluded to before, one may suspect that the samples have been subjected to some 
post-reactor contamination. To estimate this effect, we assume that a fraction £ of the 
natural abundance iV nat (see (|14D) flowed into the reactor from the outside environment, 
contributing to the observed amount of Sm. We subtract this inflow from the now observed 
isotope abundances and apply the preceding analysis to the remainder to re-calculate o^g. 
As we can see in Table [| the result is quite insensitive to ( as long as it stays around a 
few percent. This contamination range is suggested by our analysis of Gd below. 

For ( = 0.04, for example, we obtain 

0-149 = 99 ± 10 kb (30) 
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Sample 








SF84-1469 


SF84-1480 


SF84-1485 


SF84-1492 


SD.37 


N 155 (ti)(%) 


0.5006 


0.4608 


0.6065 


0.3899 


0.5915 


iVi 56 (ti)(%) 


30.03 


29.46 


30.79 


29.90 


30.37 


N 157 (h)(%) 


0.0418 


0.2641 


0.1921 


0.0505 


0.1881 


N 158 (h)(%) 


23.66 


25.07 


26.52 


23.68 


17.31 


iVieo(ti)(%) 


11.21 


12.24 


12.93 


10.95 


6.1030 


0-155 (kb) 


30.9 


16.8 


17.8 


36.7 


26.3 


0-157 (kb) 


83.3 


8.0 


14.3 


73.7 


23.3 



Table 4: Measured isotopic rations for Gd for five samples and the corresponding, calculated 
cross sections 0155 and 1x157. 

instead of (p6|), from which follows 

AE T = 2 ± 14 meV, (31) 

for either of the temperature ranges 200-400°C and 300-400°C. Although the result (|30|) 
seems worse than (|26|) , the result fl3l|) for AE T puts an even stronger limit on its time- 
variation than (B7p, because the curves corresponding to different temperatures in Fig. [I] 
happen to all intersect near AE r = 0. 

4.3 Gd data 

In analogy with our studies of the Sm isotopic compositions and estimates of neutron- 
capture cross sections tabulated in Table ||| we now study the measured Gd isotopic com- 
positions in the Oklo RZs, which are tabulated in Table |4]. To obtain this table, we used 
the fluence and effective cross-section values which are tabulated in [17] in the same way 
as we did for Sm case. By use of the formalism developed in Section |37^ we obtain the 
estimates for the neutron-capture cross sections listed in the bottom two lines of Table 

We find that the resonance-position change AE r deduced from the cross sections in 
Table f| are generally significantly different from zero in contrast to what was indicated 
by the Sm results. It is therefore natural to investigate if the effect of contamination here, 
for Gd, plays a more crucial role than for Sm, which we found was relatively insensitive 
to modest contaminations. For Gd we divide the five samples into two different classes, 
which are presented separately in Tables |5] and || The samples SF84-1469 and SF84-1492 
presented in Table |5| are characterized by much smaller abundances of 157 Gd than are the 
other three samples presented in Table [| 

We assume that for each sample the unknown contamination parameter ( is the same 
for 155 Gd and 157 Gd, and investigate if there is any £ for which we obtain the same AE r 
for both isotopes. This is based on the assumption that the mass scale M. that relates 
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SF8 


4-1469 


SF84-1492 


c 




0"157 


0155 0"l57 




(kb) 


(kb) 


(kb) (kb) 


0.000 


30.9 


83.3 


36.7 73.7 


0.001 


31.4 


102.9 


37.4 87.2 


0.002 


31.9 


134.7 


38.1 106.6 


0.003 


32.3 


195.0 


38.8 137.1 



Table 5: Calculated cross sections CJ155 and 0157 for two samples and four values of the post- 
reactor contamination £. 

AE T to Aa are the same for the two isotopes. In fact M. defined by (see (^4[) later) 



a—^=M, 32 
da 

is found to be nearly common (~ — l.lMeV) according to a phenomenological Bethe- 
Weizsacker formula, for example. Although the result might be corrected roughly by a 
factor of two if a more precise measurement is used [13], we consider the above simple 
estimate to be sufficient for our present purpose. In the same context, it seems reasonable 
to assume that M. is nearly the same also for Sm, at least sharing the same sign. 

For the sample SD.37 we have in Fig. |^ plotted AE T (Q determined in the following way. 
The neutron-capture cross sections <5- 155 and o"i 57 corresponding to the observed isotopic 
distributions with the assumed contamination removed, are calculated. The resonance- 
position changes corresponding to these cross sections are then obtained for the assumed 
reactor temperature limits from the appropriate curves in Figs. |2] and ^. 

Like in Fig. |l|, we have also two ranges. It is now the left-branch ranges that passes 





SF84-1480 


SF84-1485 


SD.37 


c 


°155 


<5"157 


<5"l55 


0157 


C"l55 °"l57 




(kb) 


(kb) 


(kb) 


(kb) 


(kb) (kb) 


0.00 


16.8 


8.0 


17.8 


14.3 


26.3 23.3 


0.01 


20.2 


11.1 


20.6 


26.4 


28.1 29.9 


0.02 


25.4 


20.6 


24.4 


236.9 


30.2 42.0 


0.03 


34.4 


227.1 


30.1 




32.7 70.9 


0.04 


53.5 




39.2 




35.6 235.3 



Table 6: Calculated cross sections CT155 and CT157 for three samples and five values of the post- 
reactor contamination £. Dash ( — ) entries are given when the calculated cross sections are 
much larger than any of those occurring in Fig. |3|. 
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0.030 0.032 0.034 0.036 0.038 0.040 
Contamination Parameter £ 

Figure 5: For the sample SD.37 we plot the resonance-position change AE r for 155 Gd and 157 Gd 
versus the contamination parameter £, the ratio of the post-reactor contamination to the natural 
abundance. The two bands show the range of AE r for the temperature range 200-400° C, with 
the lower bound of each band corresponding to T = 200°C and the upper bound to T = 400°C. 
From the location of the overlap region of the two bands we deduce that the contamination 
parameter £ is in the range 0.0355 to 0.0403 and that the resonance-position change AE r is in 
the range —26 meV to +9 meV. 

through the area covering AE T = 0. For the right-branch ranges we obtain large positive 
values for AE r . This branch is not consistent with the Sm data giving either a near- 
zero branch or a negative branch range for AE T . The requirement that the Sm and Gd 
data give the same result for the time variability of a therefore only allows the near-zero 
branches for Sm and Gd as acceptable choices for AE T . 

For the assumed temperature range 200-400°C, the bands for 155 Gd and 157 Gd corre- 
sponding to the left-branch ranges do intersect each other for ( in the range 0.036-0.040, 
yielding the range —26 meV to +9 meV for AE r , as shown in Fig. |5|. We also find that the 
band for 155 Gd is rather flat whereas that for Gd shows a rapid increase as ( changes. 
For this reason the common intersection gives AE r staying rather close to zero, even if 
the assumed equality for the two AE r 's is somewhat modified. Consider the equality is 
changed by a factor two, for example, in either direction. We find that the common range 
for ( stays near 3 ~ 4%, while that for AE r is expanded nearly twice as large, still leaving 
the conclusion of the consistency with the Sm data basically unchanged. 

Similar results are obtained for the other samples except for SF84-1492 for which no 
reasonable solution is obtained. 

In Fig. H we summarize the ranges of ( and AE r obtained for the four samples. It is 
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Figure 6: Summary of the ranges of AE r and £ estimated from Gd data from four samples, for 
the assumed temperature range T = 200-400° C. 

remarkable to find that, with the exception of SF84-1485, the obtained values of AE T are 



consistent with the Sm result (|27|). The presence of contamination is clearly indicated 
and its magnitude is also consistent with the Sm data. 

As for SF84-1485 and SF84-1492, one may suspect non-uniform contamination inside 
the reactor core or some other type of error. This suggests that more stringent limits 
on the variation of AE r would be achieved by collecting additional samples. It might 
be argued that assuming a common value of the contamination parameter ( for the two 
isotopes is not well founded. However, it is the most simple assumption, and with this 
assumption the Sm and Gd results are entirely consistent with each other, both for the 
contamination parameter ( and for the resonance-position change AE T for which the Sm 
and Gd mutually support the conclusion AE T = 0, though the conclusion is still tentative 
before the above plausibility assumptions are justified by more careful considerations, 
taking the need for more samples aside. 

We also notice that, unlike in the Sm data, the result on Gd depends on the temper- 
ature rather sensitively. If we assume twice as wide the temperature range 200 — 600°C, 
for example, the latitude of AE r in Fig. 6 extends roughly in the same rate. This tends 
to weaken the consistency check between the null results of Sm and Gd, but not to the 
extent that the non-null results are allowed to be mutually consistent. From this point of 
view, further improving the temperature estimate is strongly encouraged. 

As an additional remark about the temperature, the same argument as in Sm showing 
that temperatures higher than 1000°C tend to constrain the cross sections to too small 
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values compared with the measured values applies to Gd as well, though the result ought 
to be more complicated. 



5 Relation of the resonance-position change to the 
variation of the coupling constants 

Following Damour and Dyson we consider AE r as the change of the mass scale A4 C 
representing the difference of the expectation values of the Coulomb energy H c ; 

M c =< H c>1 -< H c > a , (33) 

where the subscripts 1 and 2 refer to the excited state of 150 Sm and the ground state of 
149 Sm plus a neutron, respectively. They estimated Ai c ps — l.lMeV, which is proportional 
to a so that [13] 

^ = — ■ (34) 
M c a 

Identifying AAi c with (|27l) we obtain 

- = (-0.2 ±0.8) x 10- 17 y _1 , (35) 
a 

which is about five times more stringent than in Ref. [13]. The improvement is mainly 
due to our choice of the separated right-branch range rather than the combined range. 

If, on the other hand, we choose (p8|) corresponding to the left-branch range, we would 
obtain 

- = (4.9 ±0.4) x KT^y" 1 , (36) 
a 

indicating an apparent evidence of the time- variability of a. 

We note that E r = 97.3meV of the resonance in question is extremely small compared 
with \M C \ ~ lMeV, most part of which is cancelled by the contribution from the strong 
interaction J\4 S defined analogously with J\4 C ; 

M c « -M s . (37) 

We assume Ai s is proportional to a™ with n likely to be 1, but allowing n ^ 1 anticipating 
some complication of nuclear physics. Combining this with (B4|) yields 



Aa s 1 Aa , , 

— « , 38 

a s n a 

if the change of the resonance energy is solely due to the change of the strong interaction 
coupling constant. In this way we expect simply the minus of (|35|) and (|36|) for the 
assumed value n = 1. These are somewhat weaker than what we might find by assuming 
M. s ps 40MeV, as suggested originally by Shlyakhter, who naively ignored the fact that 
the ground state 149 Sm has its own binding energy, which should be crucial for the near 
cancellation between the Coulomb and nuclear energies. 
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In principle, one might argue that the nuclear energy levels should be expressed in 
terms of the more fundamental constant «qcd- When applied to realistic low-energy 
phenomenology, however, it does not appear that this approach has reached the same 
level of rigor as QED. We instead expect to obtain a useful insight by expressing the 
result on the energy shift in terms of a s , which has played important roles widely in 
understanding low-energy nuclear physics. 

On the other hand, the argument can also be extended to a perhaps more realistic 
situation in which AE T is a combined effect both of the Coulomb and strong interactions. 
To show this we rewrite fl3"T| ) without omitting E r ; 

M c + M s = E T . (39) 



We also assume n — 1, for the moment. By varying a and a s in (|39|), we obtain 



— M c + ^M s = AE r . (40) 

a a s 

This is not enough to determine Aa/a or Aa s /a s separately. We need another relation 
between the two quantities. Here we simply introduce the ratio u = (Aa s /a s )/(Aa/a) 
as a convenient parameter, in terms of which we solve (EO) for Aa/a; 



= ^ . (41) 

a (1 - u)M c + uE r v 1 

Unless u is very close to 1, we find 

( l- u) ^=(iz^A* = ^, (42) 
a V u J a s M c 

hence giving d/a essentially of the same size as (|35| ) or (|36|). 
If, however, 

u - 1< — ~ O (l0~ 7 ) , (43) 



the denominator in ( |41~D is almost E r , giving 

Aa _ Aa s _ AE T 

a a s E r 

which could be as large as 0.2, resulting in d/a as large as 10 _10 y _1 . The condition 
implies, however, an unrealistic fine-tuning of parameters of any theoretical models. 



(44) 



6 Discussions 

We have re-examined in detail the derivation of the upper bound on the shift with time of 
the resonance energy of 149 Sm, and applied it to the new, carefully collected samples from 
deep- lying Oklo reactor zones. We reached the conclusion |Ai? r | <,16 meV. The surprising 
agreement to Shlyakhter's original result is probably somewhat accidental, since he in early 
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years probably had access to few uncontaminated samples and since he also assumed the 
unrealistically low temperature T = 20°C. 

We still have another choice of nonzero shift AE T ~ — lOOmeV. From the former 
more favored choice, we derived \a/a\ ~ |d s /a s | <10~ 17 y _1 , unless there is any reason 
why the equality a/a = a s /a s holds true to an extreme accuracy of 10~ 7 . The near 
equality for the time-variability of the two coupling constants is a consequence of the 
near cancellation between the two interactions resulting in an extremely small value of 
the resonance energies. 

The cross section deduced from the 15 samples analyzed in Ref. [13] is cr 149 = 75±18 kb 
(two standard deviation). We find that the corresponding critical temperature for joining 
two ranges is T crit m 900°C, which is sufficiently higher than their choice T upper = 700°C, 
implying that the two regions can be considered well separated. In this sense we may 
interpret their result in terms of the two ranges; 

-94±26meV, and 46 ± 44meV. (45) 

The separation and the sum of the half-widths are 140meV and 70meV, respectively. The 
authors seemed nevertheless to favor an even more conservative attitude to consider as if 
the two ranges were connected, the far ends of which give 

- 120meV < AE < 90meV. (46) 

The limits might be narrowed slightly to —115 meV and 68 meV, respectively, at the level 
of one standard deviation, to be compared with our analysis. Notice also that AE r = 
lies just on the edge of the two-standard deviation bound in the second of (^5|), indicating 
probably contamination of their samples. 

We also studied the isotopic composition of Gd in Oklo, to look, for the first time, 
for a possible time- variation in the position of the low-lying resonances in 155 > 157 Gd and a 
corresponding time-variation of the coupling constants. We used the same methodology 
as in the studies of Sm. In the Gd analysis it was crucial to determine the amount of 
post-reactor contamination, in contrast to the case in the Sm analysis. By assuming that 
the energy-shifts of both isotopes of Gd are nearly the same having the same sign as 
that of Sm, the results of the Gd studies confirm and reinforce the conclusion of the null 
result based on the 149 Sm data. To consolidate this tentative but reasonable conclusion, 
further analyses might be required. This should include, in addition to collecting more 
samples, re-examining the proportionality between Ai c and a as assumed in [13], because 
any change in the Coulomb energy affects the size of a nucleus hence changing the kinetic 
energy as well, resulting in rather complicated dependence of the total energy on a in 
general. This aspect might be even more relevant when we add the contribution from the 
strong interaction. 

We hope that we have shown with our detailed analysis here, having somewhat im- 
proved the result due to Damour and Dyson, that the constraints from the Oklo phe- 
nomenon should be considered to be quite realistic upper bounds on the time-variability 
of the fundamental constants. This result should provide valuable guidance in the search 
for viable unification models. 
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The result obtained here, that the upper bounds on the rate-of-change of the funda- 
mental constants are much lower than the value ~ 1CT 10 y -1 expected naturally in terms 
of the present age of the universe is likely linked to the cosmological constant problem, 
suggesting the presence of a small but nonzero A [23-27]. One of the unique features of the 
approach in this paper, in contrast to the cosmological studies of the time- variation of the 
fundamental constants, is that we probe the time ~ 2 x 10 9 y ago. This is somewhat less 
distant than many of the cosmological phenomena considered so far, which have typical 
time scales of ~ 10 10 y [3,5,8]. The different time scales of the studies may be useful in 
trying to understand the recent finding that the fine-structure constant seems to exhibit 
a complicated time dependence [28]. 
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for pointing out to us the importance of the effective cross section in the analysis of the 
Oklo phenomena. 
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Appendices 



A Numerical details 

The integral 

roo — fi? 

Jo (-B — -E r ) 2 + Itat/4 

occurring in Eq. ([8]) has to be evaluated by numerical quadrature. So that we can use 
Gauss-Laguerre integration we make the substitution x = E/T and obtain 

roo ^ roo 

I C = T dxVxTexp(-x)— J 2 , =T dxexp(-x)fi(x). (A.2) 

Jo (xT — E T ) 2 + Itot/ 4 Jo 

where 

«*> = ^V_^ +r?ot/4 - < A - 3 > 

Because /1 has a maximum at 



g L +^4^+3(rW2g 

it is necessary to divide the integration interval into two parts to achieve satisfactory 
numerical accuracy. The first part of the integration should include the maximum of the 
function f\ and thus go from to x& where 

x d = 2x max (A. 5) 

is a suitable choice. However, the numerical results are very stable to reasonable variations 
of the choice of x&. Thus we rewrite 



[ Xd I 

h = T L dxy/xfeM - x \xT-E T y+n t /A 

roo £ 

+ T dxVxTexp(-x) _ J (A.6) 



E 

{xT-^ + TljA 

So that we can use Gauss-Laguerre integration on the latter integral we substitute y = 
x — Xd and obtain 

f X d E~ 

h = Tj dxV^expi-x) (A.7) 



(xT - E r f + T t 2 ot /4 
+ T J dy^J (y + x d )Texp(-x d ) exp(-y)-^ . .. ; , , 2 p2 



Ej 



((y + x d )T-E r ) 2 + Tl t /A 

The first integral can now be calculated using Gauss-Legendre integration and the last in- 
tegral can be calculated by Gauss-Laguerre integration. Thus we obtain for the numerical 
evaluation of Eq. [A.l 



ic = E ^(^f ) + exp(-x d ) £ ^F L (yf) (A.i 

i=l i=l 
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where the superscripts G and L indicate Gauss-Legendre and Gauss-Laguerre integration, 
respectively, and where 

= TV^eM-z) (xT _ E ^ +rUi (A.9) 

and 

FHy } = T^T^f {{v + x4T E l r + rL/i (A.10) 

Note that the exponential term is retained in F G but not in F L . 
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B Details of the fluence determination in the Oklo 
reactor zones 



In this section the following notation is used: 

iV 14 3 The number of 143 Nd atoms per unit volume after the Oklo reactor shut- 
down. In the absence of contamination, in- and/or outflow of material this 
number is unchanged since the shutdown. 

./V147 The number of 147 Sm atoms per unit volume after the Oklo reactor shut- 
down. In the absence of contamination, in- and/or outflow of material this 
number is unchanged since the shutdown. 

./V235 The current number of 235 U atoms per unit volume in an Oklo reactor zone. 

./V238 The current number of 238 U atoms per unit volume in an Oklo reactor zone. 

N^ 35 The number of 235 U atoms per unit volume in an Oklo reactor zone at the 
start of the reactor about two billion years ago. This number can simply 
be calculated from data observed now at other uranium deposits where no 
nuclear reaction has occurred. 

A^2 38 The number of 238 U atoms per unit volume in an Oklo reactor zone at the 
start of the reactor about two billion years ago. This number can simply 
be calculated from data observed now. at other uranium deposits where no 
nuclear reaction has occurred. 

I a Resonance integral. 

r rz Epi-thermal index. Accounts for the contribution of the neutron spectrum 
integrated over low-lying resonances. The subscript rz indicates that de- 
pends on the shape of the neutron spectrum and other factors and is there- 
fore different in each reactor zone (rz) at Oklo. 



^143 


Cumulative fission yield of 143 Nd. 


^147 


Cumulative fission yield of 147 Sm. 


Of235 


235 U fission cross section. 


0"143 


143 Nd total neutron-absorption cross section. 


ai47 


147 Sm total neutron-absorption cross section. 


0235 


235 U total neutron-absorption cross section. 


0238 


238 U total neutron-absorption cross section. 


0°A 


Thermal neutron-capture cross section. 
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A235 235 U decay constant. 
A238 238 U decay constant. 

rz Average neutron flux in an Oklo reactor zone. 

C rz Restitution factor of 235 U in an Oklo reactor zone. It accounts for the 
restitution of 235 U through neutron capture on 238 U, /3-decay of 239 U to 
239 Np, /3-decay of 239 Np to 239 Pu, a-decay of 239 Pu to 235 U, and smaller 
contributions from other reactions. 

At rz Time duration of the neutron flux in the Oklo reactor. 

t Time elapsed since the start of the Oklo reactor, about 2 billion years. 

The fluence 4> TZ At rz in the different Oklo reactor zones may be determined from studies 
of the remnants of the reactions: 

143 Nd(n, 7 ) 144 Nd 
147 Sm(n, 7 ) 148 Sm 
235 U(n,:r)X and 
238 U(n,x)X 

The absorption cross sections a U3 and 6" 147 associated with the first two reactions may 
be expressed in the form 

a A = a A + r rz I A (B.ll) 

Specifically [29] 

<5"i43 = 335 — 100r rz 
<7 147 = 52 + 600r rz 

in units of barn. The differential equations for the time-dependence of the number of 
atoms per unit volume N A of the isotopes involved in these reactions are 



dN U3 (t) 
dt 

dN u7 (t) 
dt 

dN 235 {t) 
dt 

dN 238 (t) 
dt 



-O"l430rziVl43(£) + ^235 (O^143^f2350rz 

-vu7(f>TzNu 7 (t) + N 23 5(t)Y U7 a i235 (f) TZ 
-A 2 35^235(i) + N 235 (t)(l - C rz )a 235 r 
-X 23 8N 238 (t) + N 238 (t)a 238 (f) TZ 
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We obtain as solutions 



N- 



113 



N- 



117 



^"143 — (1 — Crz)0"235 

-^235^1470"f235 
0"147 ~~ (1 ~~ Crz)<5"235 

( N 235 



exp (-o-i4 3 rz At rz ] } (B.12) 



23S 



, ^238 . 



{exp -(1 - C rz )(3-2350rzAt r 

{exp -(1 - C rz ) CT235 rz At rz - exp (-O-i470r Z At rz J } (B.13) 
exp [-(A 23 5 - A 2 38)^] exp {- [(1 - C rz )6" 2 35 - frzss] </> rz At rz }(B.14) 



One may easily obtain the ratio N^/Nlj^ of 235 U to 238 U at Oklo by observing that 
today the ratio in uranium mines is 0.007252 and use 



, ^238 , 



exp (-A 235 t) 
exp (-A 2 38*) 



0.007252 



This leads to 



A' 



235 



N. 



238 



0.007252 exp {- [(1 - C IZ )a 235 - a 238 } rz At rz } 



(B.15) 



(B.16) 



for the ratio 235 U to 238 U observed at Oklo today. 

For each reactor zone at Oklo there are three unknowns that have to be determined: 
the restitution factor C rz , the epi-thermal index r rz , and the neutron fluence rz At rz . 
They can now be obtained as solutions to the three coupled equations ( p.l2|) , ( B.13|) , 
and ( [BlBD . Condit ions for the occurrence of fission processes are influenced by many 
factors, such as the U content, neutron-absorber elements, water content among others. 
This explains why we often find different neutron fluences in different reactors. 
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